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The  effects  of  a  two-dimensional,  cylindrically  symmetric  periodic  structure  on  the  radiation  from  a  source  lo¬ 
cated  inside  the  structure  are  considered.  The  coupling  of  the  source  radiation  to  the  cylindrical  cavity  is  ana¬ 
lyzed  classically  by  consideration  of  the  interaction  of  a  current  line  source  with  its  own  radiated  field 
supported  by  the  cavity.  The  analysis  predicts  variations  in  the  radiative  damping  rate  (inverse  lifetime)  and 
associated  shifts  in  the  oscillation  frequency  of  the  source.  It  is  found  that  significant  enhancement  and  inhi¬ 
bition  of  the  radiation  are  possible  even  when  the  source  is  not  at  the  center  of  the  structure.  Frequency  shifts 
induced  by  the  cavity  are  found  to  be  negligible  relative  to  the  size  of  the  band  gap  of  the  Bragg  structure.  The 
class  of  periodic  structures  analyzed  has  potential  application  to  planar  waveguide  devices,  such  as  concentric- 
circle  gratings,  surface-emitting  lasers,  and  cylindrically  symmetric  lasers  and  amplifiers. 


1.  INTRODUCTION 

It  is  well  understood  that  the  radiative  decay  of  excited 
atoms  through  spontaneous  emission  can  depend  strongly 
on  the  electromagnetic  modes  associated  with  the  atoms’ 
surrounding  environment.'  The  spontaneous  emission 
might  be  either  enhanced  or  inhibited,  depending  on  the 
characteristics  of  the  atom  and  the  cavity.  Many  in¬ 
vestigations  of  this  phenomenon  have  been  concerned 
with  the  effects  of  metal  and  dielectric  interfaces,  and 
metallic  waveguides  near  cutoff,  on  atomic  radiation. A 
beautiful  series  of  experiments  by  Drexhage'*  and  others 
examined  the  dependence  of  the  fluorescent  lifetime  of 
excited  molecules  on  the  separation  from  gold,  silver,  and 
copper  interfaces.  While  they  are  fundamentally  highly 
significant,  metallic  structures  are  less  important  in  prac¬ 
tice  because  of  their  inherently  lossy  nature.  As  a  result 
periodic  dielectric  structures  have  received  considerable 
attention  recently.  Most  investigations  of  these  struc¬ 
tures  have  been  concerned  with  one-dimensional  and 
three-dimensional  rectangular  periodicities  — the  one¬ 
dimensional  case  has  direct  applications  to  vertical-cavity 
surface-emitting  semiconductor  lasers;^'"’’  the  three- 
dimensional  structure  offers  the  possibility  of  a  total 
photonic  band  gap,  in  which  electromagnetic  wave  propa¬ 
gation  is  isotropically  forbidden  within  an  energy  gap  of 
nonzero  width.*  Two-dimensional  (2D)  periodic  struc¬ 
tures  have  also  received  recent  interest,*'*  In  addition  to 
their  fundamental  interest,  these  structures  may  have  ap¬ 
plications  in  planar  waveguides  and  in  lasers  and  amplifi¬ 
ers  that  operate  along  the  nonperiodic  axis  (since 
spontaneous  emission  is  inhibited  in  the  periodic  plane). 
For  model  systems  of  dimensions  higher  than  one,  previ¬ 
ous  research  focused  on  calculating  the  energy  gaps  for 
various  types  of  translationally  invariant  systems.  Of 
equal  interest,  however,  is  the  problem  of  an  elementary 
source  placed  within  a  Bragg-reflecting  structure,  a  com¬ 
plex  problem  for  practical  systems.  In  this  pajier  we  con¬ 
sider  the  possibility  of  altering  radiative  properties  in  two 
dimensions  by  using  only  a  ID  periodicity,  that  is,  along 


the  radial  direction  in  a  cylindrically  symmetric,  periodic 
structure. 

Cylindrically  symmetric,  periodic  structures  have  al¬ 
ready  received  considerable  attention  because  of  their 
potential  to  form  circularly  symmetric,  surface-emitting 
semiconductor  laser  resonators."  '*  These  lasers  employ 
a  radial  surface  corrugation  on  a  planar  waveguide.  An 
example  of  a  typical  structure  fabricated  for  use  in  a 
circular-grating  laser  is  shown  in  Fig.  1.  These  gratings 
were  defined  by  electron-beam  lithog  aphy  and  were  sub¬ 
sequently  etched  into  the  AlGaAs  laser  material  by  chemi¬ 
cally  assisted  (chlorine)  ion-beam  etching.'*  Recently  the 
first  observation  of  true  circular-mode  lasing  in  one  of 
these  lasers  was  reported.''*  For  this  case,  in  which  lasing 
occurs  along  the  radial  direction,  it  is  most  desirable  to 
enhance  spontaneous  emission,  since  it  has  been  shown 
that  one  can  decrease  the  threshold  of  a  laser  substan¬ 
tially  by  increasing  the  coupling  efficiency  of  spontaneous 
emission  into  a  lasing  mode.* 

Because  of  the  fundamentally  unique  nature  of  the 
cylindrically  symmetric  periodicity  and  the  opportunity 
for  insight  into  practical  applications  such  as  the  circular¬ 
grating  laser,  it  is  interesting  to  consider  the  cavity  ef¬ 
fects  of  such  a  structure  on  the  radiation  of  an  enclosed 
source.  At  first  thought  it  would  seem  that  the  source 
would  have  to  be  placed  at  the  center  of  the  structure  in 
order  for  it  to  see  the  periodicity  and  strongly  couple  to 
the  cavity  or  to  exhibit  significant  enhancement  or  inhibi¬ 
tion  of  radiation.  However,  for  a  cylindrical  periodic 
structure  of  infinite  radial  extent,  any  point  at  a  finite 
radius  appears  to  be  at  the  center  when  viewed  from  infin¬ 
ity;  the  radiation  from  a  source  placed  at  a  nonzero  but 
finite  radial  position  can  be  described  by  a  finite  number 
of  azimuthal  waves,  all  of  which  approach  circular  waves 
as  they  radiate  out  to  infinity.  Therefore  a  cylindrical 
periodic  structure  of  large  but  finite  radial  extent  ought  to 
couple  strongly  to  radiation  emitted  from  a  region  of 
nonzero  extent  around  the  center.  It  is  our  intent  to 
show  that  this  region  does  exist  and  can  extend  over  a 
distance  of  many  wavelengths,  indicating  that  practical 
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Fig.  1.  Scanning  electron  micrographs  of  a  typical  grating  reso¬ 
nator  for  a  concentric-circle  grating,  surface-emitting  semicon¬ 
ductor  laser. 


applications  involving  the  cylindrical  periodic  structure 
are  feasible. 

Because  the  unique  radiative  effects  associated  with 
cylindrically  symmetric,  periodic  structures  are  a  direct 
result  of  the  2D  circular  nature  of  these  structures,  we 
examine  a  strictly  2D  system  in  this  paper;  i.e.,  we  assume 
a  current  line  source  that  emits  strictly  cylindrical  waves, 
so  that  both  source  and  cavity  are  translationally  invari¬ 
ant  along  the  cylinder  axis.  This  simplification  ignores 
radiation  that  propagates  along  the  cylinder  axis.  For  a 
practical,  realizable  structure  with  cylindrical  symmetry, 
accurate  predictions  of  the  magnitudes  of  the  radiative  ef¬ 
fects  should  take  into  account  the  existence  of  radiative 
channels  associated  with  propagation  in  the  axial  direc¬ 


tion.  Thus,  while  the  theory  in  this  paper  is  not  capable 
of  quantitatively  predicting  radiative  effects  in  devices 
like  that  shown  in  Fig.  1,  it  is  likely  that  the  qualitative 
features  emerging  from  the  model  will  offer  insight  into 
the  behavior  of  such  devices.  This  is  particularly  true  for 
cylindrical  structures  that  include  a  planar  optical  wave¬ 
guide,  as  a  result  of  the  significant  coupling  between  the 
radially  propagating  bound  modes  and  a  source  located  on 
or  within  the  waveguide.'’’  In  this  context  the  theory  in 
this  paper  can  be  viewed  as  a  first  step  toward  under¬ 
standing  the  radiative  properties  of  practical  structures, 
such  as  those  involving  concentric-circle  gratings,  by  ex¬ 
amining  a  simpler  system  for  which  essentially  exact  re¬ 
sults  can  be  obtained. 

The  effect  of  a  cavity  on  the  enhancement  or  the  inhibi¬ 
tion  of  spontaneous  emission  can  be  theoretically  modeled 
in  many  ways.  A  simple  and  intuitive  model  treats  the 
source  as  a  classical  dipole  oscillator  that  is  driven  by  the 
reflection  of  its  own  radiated  field  from  the  cavity.  This 
method  was  employed  by  Kuhn  to  describe  the  effect  of  a 
metallic  mirror  on  the  radiation  from  a  nearby  excited 
molecule.'"  It  was  later  used  quite  successfully  by 
Chance  et  al.''  to  model  quantitatively  the  experimental 
data  of  Drexhage.  We  adapt  this  classical  method  to  ana¬ 
lyze  the  effects  of  the  cylindrical  periodic  structure  on  the 
radiation  of  an  enclosed  source  for  several  reasons.  First, 
the  method  allows  us  directly  to  analyze  radiative  effects 
as  a  function  of  source  position  in  the  cavity  and  hence  to 
address  questions  regarding  the  size  of  the  central  region 
for  which  significant  effects  can  occur.  Second,  it  is  not 
clear  how  a  band-theory  description  of  electromagnetic 
wave  propagation  can  be  applied  to  the  cylindrical  case,  as 
is  often  done  in  the  case  of  the  rectangular  periodicity. 
In  particular,  many  concepts  in  band  theory,  such  as  re¬ 
ciprocal  space  and  the  Brillouin  zone,  rely  on  the  transla¬ 
tional  invariance  inherent  in  rectangular  periodic  systems. 
In  the  cylindrical  system,  if  one  moves  an  integer  number 
of  periods  along  the  radius,  the  surrounding  structure 
changes,  and  hence  the  system  is  not  translationally 
invariant. 

The  organization  of  this  paper  is  as  follows.  I;;- 
Section  2  we  develop  the  theoretical  formulation  for  calcu¬ 
lation  of  the  enhancement  and  inhibition  of  cylindrical- 
wave  radiation  from  a  current  line  source  in  terms  of  the 
radiative  damping  rate  (inverse  lifetime)  and  the  associ¬ 
ated  frequency  shifts.  The  formulation  uses  cylindrical 
coupled-wave  theory  to  calculate  the  reflection  of  the 
radiated  field  back  onto  the  source."  '"  In  Section  3  the 
theory  is  used  to  calculate  numerically  radiation  enhance¬ 
ment  and  inhibition  and  frequency  shifts  for  representa¬ 
tive  structures.  Finally,  in  Section  4  conclusions  are 
drawn  from  the  results  of  this  investigation. 

2.  THEORETICAL  FORMULATION 

In  this  analysis  we  treat  the  source  of  radiation  as  a  classi¬ 
cal  current  line  source.  We  model  the  effects  of  the  cylin¬ 
drical  periodic  structure  on  the  source,  following  Kuhn’s 
method.'®  In  Subsection  2.A  we  describe  the  equation  of 
motion  of  the  current  line  source,  and  from  that  equation 
we  derive  the  changes  in  radiative  damping  rate  (inverse 
lifetime)  and  frequency  from  the  free-space  quantities 
that  are  due  to  the  source’s  own  field  reflected  from  the 
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cavity.  In  order  to  calculate  the  reflected  field,  we  must 
first  calculate  the  field  emitted  by  the  source;  this  field  is 
derived  in  Subsection  2.B.  In  Subsection  2.C  we  derive 
the  reflected  field  at  the  source  position  by  using  cylindri¬ 
cal  coupled-wave  theory.  Finally,  in  Subsection  2.D,  com¬ 
plete  expressions  for  the  radiative  damping  rate  and  the 
frequency  shift  are  derived. 


A.  Source  Equation  of  Motion 

In  Kuhn’s  analysis  the  source  is  assumed  to  be  a  harmoni¬ 
cally  bound  charge  that  represents  a  point  dipole.  In  the 
present  2D  analysis  the  reduction  in  the  order  of  the  prob¬ 
lem  dimensionality  must  be  accompanied  by  a  redefinition 
of  characteristic  source  parameters  in  terms  of  a  unit 
length.  We  assume  that  the  current  line  source  behaves 
as  an  ensemble  of  harmonically  bound  charge  elements, 
described  by  the  charge  per  unit  length,  q,  such  that  this 
ensemble  obeys  the  equation  of  motion 

p  -I-  aia^p  =  (q‘^/m)ER  -  bap,  (1) 

where  p  is  the  total  dipole  moment  per  unit  length,  tuo  is 
the  free-space  oscillation  frequency  in  the  absence  of  all 
damping,  m  is  the  effective  mass  per  unit  length,  Er  is  the 
reflected  field  at  the  source  position,  and  6o  is  the  damp¬ 
ing  constant  (inverse  lifetime)  in  the  absence  of  the  cavity. 
The  dipole  moment  per  unit  length,  p,  and  the  reflected 
field,  Er,  oscillate  at  the  same  (complex)  frequency: 

p  =  Po  exp[-(i6j  +  b/2)t\ ,  (2a) 

Er  =  Eo  exp[-(i<a  +  b/2)t],  (2b) 

where  w  and  b  are  the  frequency  and  the  damping  rate  in 
the  presence  of  the  cavity,  respectively.  The  problem  is 
then  one  of  a  driven  harmonic  oscillator,  where  the  exter¬ 
nal  force  is  proportional  to  the  reflected  radiation  field  of 
the  line  source. 

Using  the  forms  given  by  Eqs.  (2)  in  Eq.  (1)  and  equating 
real  and  imaginary  parts  independently,  we  arrive  at  the 
expressions 


-  =  1  -f 

■■■■  ■  Im(Eo), 

(3) 

^0 

/TZWPoOo 

2  2 

0)  -  (jjf)  =  -r  - 

■^--^Bb(Eo). 

(4) 

Equation  (3)  gives  the  normalized  (to  the  free-space  value) 
radiative  damping  rate.  Equation  (4),  which  gives  the  fre¬ 
quency  shift  caused  by  the  presence  of  the  cavity,  can  be 
simplified  by  recognizing  that  A<u  «  <o,  <wo,  where  = 
<0-0)!,.  Using  this  assumption,  which  will  be  justified 
below,  we  can  write  Eq.  (4)  as 


Ao)  = 


8a)o 


66o 

4^0 


2 


2fTl(Oo  Po 


Re(Eo)  • 


(5) 


In  order  to  calculate  the  radiative  damping  rate  and  fre¬ 
quency  shift  determined  by  Eqs.  (3)  and  (5),  respectively, 
we  need  to  know  the  reflected  field  strength  at  the  source, 
Eo-  It  is  convcriept  to  factor  this  quantity,  such  that 


Eo  —  ErF.  (6) 

The  quantity  Es,  which  denotes  the  source  field  strength, 
will  be  calculated  in  Subsection  2.B  of  this  section.  The 


quantity  F  is  an  effective  reflection  coefficient  that  de¬ 
pends  on  the  source  position  and  the  characteristics  of 
the  cylindrical  periodic  structure.  In  Subsection  2.C  we 
show  how  this  quantity  can  be  computed  with  coupled- 
wave  theory. 


B.  Source  Field 

Reducing  the  problem  to  one  of  strictly  2D  character,  we 
take  the  source  to  be  a  current  line  source  oriented  along 
the  cylinder  axis  of  the  cylindrical  periodic  structure. 
The  geometry  of  the  problem  is  shown  in  Fig.  2.  The 
source  is  described  by  the  current  density 

J(r,(j>)  =  -ioJoPoStr  -  ro)z 

=  -!ajoPo(Vr)5(r  -  ro)S{<l>  -  4>o)^,  <7) 

where  ro  =  (ro,4>o)  is  the  location  of  the  source  and  5(r)  is 
the  2D  Dirac  delta  function. 

Because  of  the  simple  geometry  of  the  source  and  the 
periodic  structure,  we  need  to  consider  only  a  single  vec¬ 
tor  component  of  the  electric  field,  E  =  zE,(r,<t>),  which  is 
a  solution  to  the  wave  equation 

V  X  V  X  E(r,<^)  -  ^ioe«ao^E(r,<l))  =  .  (8) 

where  an  exp(-ia(ot)  harmonic  time  dependence  is  as¬ 
sumed  and  e  is  the  constant  permittivity  in  the  source  re¬ 
gion  (r  <  rj).  From  the  fact  that  the  electric  field  has 
only  a  z  component,  Eq.  (8)  is  simplified  to  a  2D  Helmholtz 
equation  with  a  point-source  driving  term: 

(V^  -t-  ko^)EAr,<l>)  =  -poi^oPobir  -  ro),  (9) 

where  ho  =  (/^os)'  ^a>o-  The  solution  to  such  an  equation 
is  well  known  to  be 

£’i(r,d>)  =  £.5Hi)‘’(*o|r  -  Tol),  (10) 

where 


Es  ^ 


ipooio 

4eoC^ 


(11) 


and  Ho'  is  a  Hankel  function  of  the  first  kind,  order  zero. 


Fig.  2.  Illustration  of  the  geometry  of  the  current  line  source 
and  the  cylindrical  periodic  structure.  Both  source  and  struc¬ 
ture  extend  uniformly  to  z  =  ±<». 
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Although  Eq.  (10)  accurately  describes  the  field  emitted 
by  the  line  source,  it  is  not  the  most  convenient  form.  By 
expanding  the  4)  dependence  of  the  field  in  a  Fourier  se¬ 
ries,  one  can  express  the  source  field  as 


E.^(r,<i))  =  Es  2  exp|iOT(d>  -  <^o)| 

in  —  -  X 

^  ^J„,ikor)Hl]'ik„ru) 
[</m(Aorii)//;„"(/tor) 

r  <  r.) 

(12) 

r  >  r„ 

where  J„,  is  a  Bessel  function  of  order  m. 
source  field  can  be  compactly  written  as 

For  r  >  r„,  the 

E/ir,(t>)  =  Es  S  C„,H'miknr)expiim(t>) 

nj  -  X 

r  >  r»,  (13) 

with 

Cm  s  exp(-im<t>»)J„Akaro) . 

(14) 

Equations  (13)  and  (14)  clearly  demonstrate  the  nature  of 
the  source  field  emitted  by  the  line  source.  For  r  >  n,, 
the  field  is  composed  of  a  sum  of  individual,  outward¬ 
going  azimuthal  waves  of  order  m,  whose  strengths  are 
given  by  the  amplitudes  C,„ .  Even  though  the  sum  is  infi¬ 
nite,  for  a  finite  source  position  r,,  the  amplitudes  C„  are 
appreciable  only  for  a  finite  number  of  waves;  therefore 
the  infinite  sum  can  be  approximated  well  by  a  finite 
number  of  terms.  Note  that  the  source  field  is  composed 
of  standing  waves  for  r  <  ro  and  hence  is  well  behaved  at 
the  origin  for  a  finite,  but  nonzero,  source  position. 

C.  Reflected  Field 

In  this  subsection  we  calculate  the  effective  reflection  co¬ 
efficient  r  defined  in  Subsection  2.A.  The  basic  idea  is  to 
express  the  source  field  as  a  sum  of  individual  azimuthal 
wave  components  (as  in  Subsection  2.B),  to  calculate  the 
strengths  of  those  components  at  the  inner  radius  of 
the  cylindrical  periodic  structure  (r  =  ri).  to  compute  the 
Bragg  reflection  of  each  azimuthal  component  separately 
by  using  cylindrical  coupled-wave  theory,  and  then  to  sum 
the  reflected  components  back  at  the  source  to  obtain  the 
total  reflected  field. 

In  the  source  region,  defined  such  that  r  <  r,,  the  total 
electric  field  E,^  can  be  expressed  as  the  sum  of  the  source 
field  E/'^  and  a  reflected  field  E^'^,  consisting  of  outward- 
and  inward-going  cylindrical  waves  generated  by  reflec¬ 
tions  off  the  cylindrical  Bragg  structure: 

E/==E/  +  E/,  (15) 

with  E/  given  by  Eq.  (12)  and 

£/ =  2  [a„I/!^'(knr)  +  b„If!^'(ki,r)]exp(im4>),  (16) 

m  =  —  * 

where  a„  and  b„  are  (constant)  outward-  and  inward-going 
wave  amplitudes.  The  amplitudes  a„  and  can  be  de¬ 
termined  from  the  boundary  conditions  on  the  source  re¬ 
gion.  At  the  center,  r  =  0,  the  sum  of  all  outward-going 
waves  of  order  m  must  be  equal  to  the  sum  of  all  inward- 
going  waves  of  order  m.  This  intuitive  boundary  condi¬ 
tion,  indicating  perfect  reflection  at  the  center,  is 
mathematically  necessary  to  prevent  the  total  field  from 
diverging  at  the  center  owing  to  the  singular  nature  of 
cylindrical  Neumann  functions.'®  Since  this  boundary 
condition  is  already  built  into  the  source  term  of  the  total 


field,  application  to  Eq.  (16)  requires  that 

a„,=b,„.  (17) 

The  other  boundary  condition  is  applied  at  the  inner  radius 
of  the  cylindrical  Bragg  structure,  r  =  r,.  It  requires 
that  the  sum  of  all  inward-going  waves  of  order  m  be  equal 
to  the  sum  of  all  outward-going  waves  of  order  m  multi¬ 
plied  by  the  amplitude  reflection  coefficient,  p„,  of  the 
Bragg  structure  for  a  wave  of  order  m  : 

b„H'„i\karO  =  +  E.sC„, //„!  (^,jr, )| . 

(18) 


The  reflection  coefficient  p,„  is  calculated  by  using 
scalar,  cylindrical  coupled-wave  theory."*  This  develop¬ 
ment  uses  a  field  expansion  of  the  form  of  Eq.  (16),  where 
the  outward-  and  inward-going  wave  amplitudes  become 
functions  of  the  radial  coordinate  to  allow  for  the  coupling 
induced  by  the  periodic  structure  in  the  region  r,  < 
r  <  r-i.  The  coupled-wave  equations  for  these  outward- 
and  inward-going  field  amplitudes  in  the  structure,  A„,(r) 
respectively,  are  given  by 

X  \A„(r)H'„]'{kor}  +  ,  (19a) 

X  lA„(r)H!^'lkor)  +  .  (19b) 


and  B„(r), 
dA„ 
dr 


dB„ 

dr 


These  equations  are  solved  numerically  in  the  periodic  re¬ 
gion  with  the  boundary  conditions  Am  (ri)  =  1  and  Bm'Ca)  = 
0,  since  no  inward-going  waves  are  incident  upon  the 
structure  from  infinity.  In  these  equations 


Kir)  =  2k  cosi 


(20) 


where  A  is  the  period  of  the  Bragg  structure,  fl  is  the 
phase  of  the  periodicity  with  respect  to  the  center,  and  the 
coupling  constant  k  is  given  by 

K  =  -^A(n‘).  (21) 

4n 


In  this  expression  n  is  the  bulk  index  of  refraction  of  the 
source  region  (r  <  r, )  and  A(n^)  is  the  index  modulation  in 
the  periodic  region  (ri  <  r  <  r^).  The  reflection  coeffi¬ 
cient  is  then  given  by 


Pm 


B„(ri) 

A„(r,) 


(22) 


Equations  (17)  and  (18)  can  be  solved  simultaneously  to 
yield  the  constant  amplitudes  a„  and  b„.  These  can  then 
be  inserted  into  Eq.  (16),  giving  the  reflected  field  every¬ 
where  in  the  source  region.  Evaluating  this  expression  at 
the  source  position  and  recognizing  that  £o  =  E.-*(ro,<to)  = 
FEs.  we  obtain 


r=S 


2p„H':i\k,r,) 


H'^XkorO  -  p„H'J,’(kori) 


[«/m(^oro)l^ 


(23) 


for  the  effective  reflection  coefficient  of  the  radiated  field 
at  the  source  position.  For  a  given  ro,  the  Bessel  function 
J„  will  eventually  tend  to  zero  for  large  enough  m,  so  that 
the  infinite  sum  in  Eq.  (23)  can  be  well  approximated  by  a 
finite  sum. 
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D.  Radiative  Damping  Rate  and  Frequency  Shift 

From  Eq.  (3)  it  is  apparent  that,  to  calculate  the  actual 
damping  rate  b,  we  need  to  know  the  free-space  damping 
rate  6o  for  the  current  line  source.  Because  the  line 
source  in  free  space  is  radiating  energy,  its  oscillation 
must  be  damped.  One  can  calculate  the  rate  of  damping 
approximately  by  using  a  classical  radiation  reaction  argu¬ 
ment.'^"  Starting  with  the  electric  field  given  in  Eq.  (10), 
the  far-field  Poynting  vector  for  the  line  source  in  free 
space  can  be  calculated  and  then  integrated  over  the  azi¬ 
muthal  angle  to  yield  the  power  per  unit  length  emitted  by 
the  source: 


P  =  —\Es\^ 

Moa>o 


(24) 


By  application  of  the  usual  radiation-reaction  argument, 
which  is  predicated  on  balancing  the  rate  of  enei^y  lost  to 
radiation  with  the  rate  of  work  done  on  the  source  by  a 
reaction  force,  an  approximate  damping  term  is  incorpo¬ 
rated  into  the  source  equation  of  motion.  This  develop¬ 
ment  yields  an  explicit  expression  for  the  free-space 
damping  rate: 


bo 


4enmc^ 


(25) 


It  is  important  to  note  that  this  damping  rate  is  directly 
proportional  to  the  frequency  u>o,  whereas  the  damping 
rate  for  a  point  dipole  in  free  space  is  proportional  to  the 
square  of  the  frequency.  This  difference  is  representative 
of  a  fundamental  scaling  because  of  the  reduced  dimen¬ 
sionality  of  the  problem. 

Substituting  Eqs.  (25),  (11),  and  (6)  into  Eq.  (3),  we  ar¬ 
rive  at  the  simple  expression  for  the  radiative  damping 
rate  in  the  presence  of  the  cavity 

b/bo  =  1  -t-  Re(r) ,  (26) 


where  the  approximation  w  ==  a»o  is  used.  Because  ReiF) 
spans  the  range  (-l,=o),  the  normalized  damping  rate  b/bo 
will  span  the  physically  meaningful  range  (0, »).  Equa¬ 
tion  (26)  is  the  expression  that  we  use  in  Section  3  to 
calculate  radiative  damping  rates  numerically.  This  ex¬ 
pression  clearly  shows  that  the  damping  rate  is  propor¬ 
tional  to  reflected  field  strength.  Therefore,  when 
b/ba  0,  or  when  there  is  a  field  node  at  the  source  posi¬ 
tion,  the  radiative  lifetime  becomes  large,  indicating  inhi¬ 
bition  of  the  radiation.  As  b/bo  — »  or  when  there  is  an 
antinode,  the  radiative  lifetime  becomes  small,  indicating 
enhancement. 

The  frequency  shift  can  be  calculated  from  Eq.  (5)  by 
the  use  of  Eqs.  (26),  (25),  (11),  and  (6),  yielding 

L  2  L 

»  ^{[Re(r)r  -  1}  +  TT  Im(r) .  (27) 

Note  that  the  classical  frequency  shift  exactly  agrees  with 
the  quantum-mechanically  calculated  level  shift  for  the 
first  excited  state  when  the  source  is  described  as  a  har¬ 
monic  oscillator.^'  For  a  more  realistic  model  of  an 
atomic  source,  one  needs  a  thorough  quantum-mechanical 
calculation  to  predict  the  level  shift  accurately.  To 
examine  the  magnitude  of  the  frequency  shift,  consider 
a  harmonically  oscillating  line  source  emitting  at  a 
free-space  wavelength  of  Ao  =  0.5  /xm.  In  this  case 


6o  =  0.9  X  10“  Hz  and  wo  =  0.4  x  1.0'“  Hz.  Therefore 
6o  <s:  u>o,  and  relation  (27)  is  well  approximated  by 

=  (5o/2)Im(r) .  (28) 

Thus,  for  most  of  the  range  of  F  Aw  ~  bo,  indicating  that 
Aw  «  Wo,  as  we  have  assumed  above.  The  possibility  of 
large  frequency  shifts  is  addressed  more  carefully  in 
Section  3  when  we  consider  numerical  results. 

3.  NUMERICAL  RESULTS  AND  DISCUSSION 

In  this  section  the  capability  of  a  cylindrical  Bragg  struc¬ 
ture  to  enhance  or  to  inhibit  source  radiation  is  investi¬ 
gated  numerically  for  some  representative  cases.  We  are 
concerned  mostly  with  the  radiative  damping  rate,  b,  as  a 
function  of  various  parameters,  such  as  source  position 
and  the  modulation,  the  placement,  and  the  size  of  the 
Bragg  structure.  Most  of  the  cases  considered  below 
involve  structures  with  only  10  Bragg  periods,  which  is 
perhaps  the  lower  limit  that  one  might  employ  in  any  prac¬ 
tical  structure,  since  significant  reflectivity  from  a  10- 
period  grating  requires  large  modulation,  or  coupling  con¬ 
stant  K.  While  this  small  number  of  periods  was  chosen 
mostly  for  convenience  and  for  best  illustration  of  the  ra¬ 
diative  effects,  it  is  important  to  note  that  there  is  nothing 
in  the  theoretical  formulation  that  limits  the  amount  of 
modulation  for  which  the  theory  is  accurate.  Even 
though  we  employ  coupled-wave  theory,  the  theory  is  not 
perturbative  (as  in  the  case  of  an  actual  waveguide  grat¬ 
ing),  and  we  have  made  no  slowly  varying  envelope  ap¬ 
proximation.  Also,  most  of  the  cases  below  place  the 
beginning  of  the  Bragg  structure  at  least  10  periods  from 
the  center.  This  choice  allows  us  to  examine  a  reasonably 
large  range  of  source  positions,  since  the  developed  theory 
inherently  requires  the  source  to  be  inside  the  inner  ra¬ 
dius  of  the  Bragg  region  (ro  <  rj).  However,  a  more  gen¬ 
eral  treatment  could  be  extended  to  include  the  possibility 
of  source  placement  within  the  Bragg  region. 

We  first  examine  the  dependence  of  the  radiative  damp¬ 
ing  rate  on  the  phase  of  the  cosinusoidal  periodicity.  We 
assume  that  the  source  is  exactly  tuned  to  the  Bragg  fre¬ 
quency,  or  woO  =  Wfl,  where  wg  =  ttc/A.  The  phase  fi  is 
defined  by  Eq.  (20)  and  is  shown  schematically  in  Fig.  2. 
This  dependence  is  plotted  in  Fig.  3  for  the  source  posi¬ 
tioned  exactly  at  the  center  (tq  =  0).  The  Bragg  region 
extends  10  periods  along  the  radial  direction  and  consists 
of  three  different  modulation  strengths,  denoted  by  the 
coupling-constant-grating-length  products  kL  of  1,  2,  and 

4.  Over  most  of  the  full  range  of  the  grating  phase  fl  the 
source  experiences  inhibited  radiation,  since  b/bo  <  1;  the 
stronger  the  modulation,  the  greater  the  degree  of  inhibi¬ 
tion.  Near  fl  =  ■n/2,  or  for  a  sinusoidal  grating,  the 
source  experiences  strongly  enhanced  emission.  Note 
that  the  enhancement  resonance  becomes  significantly 
sharper  for  stronger  modulations.  Also  shown  in  Fig.  3  is 
the  case  of  a  100-period  Bragg  region  with  modulation 
strength  kL  =  2.  The  behavior  is  similar  to  that  of  the 
10-period  case  except  for  a  slight  shift  in  the  enhancement 
resonance. 

In  Fig.  4  we  examine  the  dependence  of  radiative  damp¬ 
ing  rate  on  the  source  position.  To  examine  the  size  of 
the  source  region  over  which  significant  enhancement  can 
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GRATING  PHASE  ( O ' "  ) 

Fig.  3.  Normalized  radiative  damping  rate  versus  grating  phase 
for  a  full  range  of  phases  (0  — •  2ir)  and  source  placement  at  the 
center. 


occur,  we  fix  ihe  grating  phase  at  the  peak  of  the  enhance¬ 
ment  resonance  in  Fig.  3  (for  ro  =  0).  This  phase  is 
denoted  ft  =  ftma«>  where  ft„,u,,  is  a  function  of  the  modu¬ 
lation  strength  kL.  Figure  4(a)  shows  the  radiative 
damping  rate  for  all  possible  source  positions  inside  the 
inner  radius  of  the  Bragg  structure  (ro  <  ri).  It  is  appar¬ 
ent  that,  no  matter  where  the  source  is  placed,  the  source 
experiences  significant  enhancement  for  this  particular 
choice  of  grating  phase.  The  largest  enhancement  occurs 
at  the  center,  since  the  source  emits  a  perfectly  cylindri- 
cally  symmetric  field  from  that  position,  and  then  the  en¬ 
hancement  decreases  as  the  source  is  moved  away  from 
the  center.  Although  the  strong  modulation  case  (kL  = 
4)  experiences  the  greatest  degree  of  enhancement,  the 
enhancement  oscillates  dramatically  with  source  position 
because  of  the  sharpness  of  its  resonance. 

In  Figs.  4(b)  and  4(c)  we  examine  the  size  of  the  source 
region  over  which  significant  inhibition  can  occur.  In 
Fig.  4(b)  we  fix  the  grating  phase  at  the  point  of  minimum 
radiative  damping  (maximum  inhibition)  in  Fig.  3.  This 
phase  is  denoted  ft  =  ft  min-  The  radiative  damping  rate 
is  plotted  for  all  possible  source  positions  inside  the  inner 
radius  of  the  Bragg  structure.  The  plot  shows  that  the 
source  experiences  significant  inhibition  for  most  source 
positions.  In  fact  there  is  a  region  with  a  radius  of  ap¬ 
proximately  two  full  grating  periods  for  which  the  source 
experiences  essentially  the  same  inhibition  as  it  does 
when  placed  exactly  at  the  center.  As  in  Fig.  4(a), 
stronger  modulation  produces  a  greater  degree  of  inhibi¬ 
tion  and  greater  oscillation  of  inhibition  with  source  posi¬ 
tion.  This  oscillation  is  not  absolutely  necessary,  though. 
In  Fig.  4(c)  we  fix  the  grating  phase  at  ft  =  Stt/S,  which  is 
near  the  maximum  inhibition  phase  but  not  exactly  equal 
to  it,  and  we  perform  the  same  calculation  as  that  in 
Fig.  4(b).  We  find  that  the  degree  of  inhibition  is  in  gen¬ 
eral  greater  over  the  range  of  source  positions  than  that 
for  ft  =  ft  rain,  and,  more  important,  the  inhibition  does 
not  oscillate  nearly  as  much  with  source  position. 

Also  shown  in  Fig.  4  is  the  case  of  a  100-period  Bragg 
region  for  kL  —  2  (dashed  curves).  The  behavior  is  simi¬ 
lar  to  that  of  the  10-period  case  except  that  the  enhance- 


(a) 


(b) 


(c) 

Fig.  4.  Normalized  radiative  damping  rate  versus  source  posi¬ 
tion  for  various  index  modulation  strengths:  (a)  grating  phase  ft 
chosen  to  give  maximum  damping  rate  (maximized  enhance¬ 
ment)  in  Fig.  3,  (b)  grating  phase  ft  chosen  to  give  minimum 
damping  rate  (maximized  inhibition)  in  Fig.  3,  (c)  grating  phase 
fi  -  37r/2. 
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SOURCE  POSITION  (fo/A) 

Fig.  5.  Normalized  radiative  damping  rate  versus  source  posi¬ 
tion  for  various  inner  radii  of  a  10-period  Bragg  region. 

merit  and  the  inhibition  remain  large  for  a  greater  range 
of  source  positions  than  in  the  10-period  case.  The  rea¬ 
son  for  this  slight  difference  can  be  explained  as  follows. 
The  field  emitted  from  a  given  source  position  can  be  well 
described  by  a  finite  number  of  azimuthal  waves,  as  im¬ 
plied  in  Eqs,  (13)  and  (14),  Over  the  10-period  Bragg  re¬ 
gion  the  highest-order  waves  are  still  somewhat  not 
cylindrically  symmetric  and  hence  are  not  strongly  re¬ 
flected  by  the  structure.  However,  because  the  100- 
period  Bragg  region  covers  a  much  larger  radial  extent, 
the  highest-order  waves  can  propagate  further  out  radially 
and  more  closely  approach  cylindrically  symmetric  waves. 
Hence  these  waves  are  more  strongly  reflected  in  the  100- 
period  case,  even  though  the  coupling-constant-grating- 
length  product  kL  is  the  same  in  both  cases. 

This  phenomenon  is  even  more  apparent  in  Fig.  5,  where 
the  inner  radius  of  the  Bragg  structure  (rO  is  varied  for  a 
fixed  grating  length.  In  Fig.  5  the  radiative  damping  rate 
is  plotted  as  a  function  of  source  position  for  three  differ¬ 
ent  inner  radii  of  a  10-period  Bragg  structure.  In  this 
plot  the  grating  phase  is  chosen  to  give  maximum  inhibi¬ 
tion  (fl  =  Stt/2).  The  farther  from  the  center  the  Bragg 
structure  is  placed,  the  larger  the  radial  extent  of  the  re¬ 
gion  over  which  the  source  experiences  maximum  inhibi¬ 
tion.  Moving  the  inner  radius  of  the  Bragg  structure 
from  10  to  40  periods  from  the  center  almost  doubles  the 
radius  of  the  maximum  inhibition  region  from  2  to  4  grat¬ 
ing  periods.  The  reason  for  this  behavior  was  explained 
above:  the  farther  from  the  center  the  Bragg  reflection 
is  occurring,  the  more  of  a  chance  the  source  waves  have 
to  approach  their  asymptotic,  cylindrically  symmetric 
form.  Hence,  for  either  enhancement  or  inhibition  of  the 
source  radiation,  it  is  desirable  to  push  Bragg  reflections 
as  far  from  the  center  as  possible,  either  by  extending  the 
inner  radius  of  a  strongly  modulated,  short  grating  or  by 
using  a  more  weakly  modulated,  long  grating.  In  a  prac¬ 
tical  .structure  the  maximum  radial  extent  of  the  Bragg 
cavity  would  be  limited  by  propagation  loss,  which  we  do 
not  consider  here. 

In  Figs.  6  and  7  we  turn  to  the  issue  of  frequency  shifts 
caused  by  the  presence  of  the  cylindrical  Bragg  cavity.  A 
natural  question  to  consider  when  using  a  Bragg  cavity. 
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whose  reflectivity  properties  are  strongly  dependent  on 
frequency,  is  whether  it  is  possible  for  the  cavity  to  shift 
the  frequency  of  the  line  source  enough  to  push  it  outside 
the  band  gap  of  the  cavity.  If  such  a  shift  were  possible, 
then  the  above  numerical  analysis,  which  assumes  that  the 
source  is  always  tuned  to  the  Bragg  frequency,  would  be 
incorrect.  The  correct  analysis,  accounting  explicitly  for 
the  interdependence  of  the  lifetime  and  frequency  varia¬ 
tions,  would  be  more  complicated.  Fortunately,  as  we 
show  below,  it  is  an  excellent  approximation  to  assume 
that  the  cavity-induced  frequency  shift  is  negligible  as  far 
as  Bragg-ref  lectivity  effects  are  concerned. 

We  first  consider  the  question  of  how  rapidly  Bragg  re¬ 
flectivity  effects  vary  with  frequency  detuning.  Thus  far 
in  the  analysis  we  have  assumed  that  the  .source  is  exactly 
tuned  to  the  Bragg  frequency,  or  =  (Ug.  This  choice 
is  made  in  the  calculation  of  the  field  reflection  coeffi¬ 
cient  from  the  Bragg  structure  with  the  coupled-wave 
equations  (19).  In  Fig.  6  we  plot  the  radiative  damping 


DETUN)NG  (o>-(i)g)/c  (cm  ') 

Fig.  6.  Normalized  radiative  damping  rate  versus  detuning  of 
the  source  frequency  from  the  Bragg  frequency  defined  by  the 
Bragg  structure.  The  line  source  is  placed  at  the  center,  and  the 
grating  phase  is  chosen  to  give  inhibition  at  the  Bragg  frequency. 
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GRATING  PHASE  (G/Jt) 

Fig.  7.  Cavity-induced  frequency  shift  of  the  line  source  versus 
grating  phase  for  a  full  range  of  phases  (0  — *  27r)  and  source 
placement  at  the  center. 


398  J.  Opt.  Soc.  Am.  B/ Vo).  10,  No.  2/F’ebruary  1993 


Erdogan  et  al. 


rate  as  a  function  of  frequency  detuning  from  the  Bragg 
frequency.  We  assume  a  Bragg  structure  designed  for  a 
source  emitting  at  Ao/n  =  0.5  /tm,  so  that  A  =  0.25  /xm, 
and  0)8  is  fixed;  we  vary  the  source  frequency  <o  =  won. 
The  source  is  placed  at  the  center  (r,)  =  0),  and  the  grat¬ 
ing  phase  fl  is  chosen  to  give  maximum  inhibition  at  the 
Bragg  frequency.  The  plot  shows  expected  oscillations 
in  the  damping  rate,  with  the  maximum  inhibition  occur¬ 
ring  at  the  Bragg  frequency.  Note  that  the  cavity  reso¬ 
nances  that  produce  local  maxima  in  inhibition  are  more 
closely  spaced  for  the  case  of  weaker  modulation  (kL  =  1), 
since  the  radial  extent  of  the  cavity  is  effectively  larger. 
It  is  important  to  note  the  frequency  scale  in  this  plot; 
the  radiative  damping  rate  does  not  change  significantly 
over  a  >1000-cm  '  range  of  frequencies  around  the  Bragg 
frequency. 

In  Fig.  7  we  examine  typical  frequency  shifts  that  are 
induced  by  the  presence  of  the  cavity.  The  plot  is  similar 
to  Fig.  3,  except  that  we  show  variation  in  frequency  shift, 
Aw/c,  instead  of  damping  rate,  as  a  function  of  grating 
phase  n.  The  source  is  again  placed  exactly  at  the  center, 
so  that  the  plot  gives  a  good  indication  of  the  full  range  of 
possible  frequency  shifts  for  any  source  position.  The 
frequency  shifts  are  generally  tiny,  except  in  the  vicinity 
of  the  enhancement  resonance  around  ft  =  Tr/2.  Even  at 
resonance  the  maximum  frequency  shift  never  exceeds 
i0.05  cm  ',  which  is  considerably  less  than  the  1000-cm  ' 
range  over  which  Bragg  detuning  can  be  considered  negli¬ 
gible.  Therefore  the  assumption  that  the  source  always 
remains  tuned  to  the  Bragg  frequency  is  a  good  one. 

4.  CONCLUSIONS 

It  has  been  shown  that  it  is  indeed  possible  either  to  en¬ 
hance  Or  to  inhibit  the  emission  of  radiation  from  a  cur¬ 
rent  line  source  located  inside  a  cylindrically  symmetric 
periodic  structure.  Whether  the  source  experiences  en¬ 
hancement  or  inhibition  depends  on  the  phase  of  the  peri¬ 
odic  grating  with  respect  to  the  center  of  the  cylindrical 
system.  When  the  line  source  is  placed  exactly  at  the  cen¬ 
ter  of  the  system,  it  tends  to  experience  inhibition  for 
most  of  the  range  of  grating  phases;  only  a  very  narrow 
phase  range  leads  to  enhancement.  This  result  indicates 
that  fabrication  of  a  practical  device  that  is  designed  to 
inhibit  radiation  will  probably  be  easier  than  fabrication 
of  an  enhancement  device  because  of  tighter  tolerances  in 
the  required  grating  phase. 

Perhaps  even  more  interesting  is  the  result  that  signifi¬ 
cant  coupling  of  source  radiation  to  the  cylindrical  Bragg 
cavity  occurs  not  only  for  placement  of  a  line  source  at  the 
center  but  also  for  a  substantial  region  of  source  positions 
in  the  vicinity  of  the  center.  This  result  is  important  for 
practical  reasons,  since  it  suggests  that  spontaneous  emis¬ 
sion  from  sources  located  anywhere  in  an  appreciable 
cross-sectional  area  of  a  cylindrically  symmetric  device 
can  be  enhanced  or  inhibited  by  cylindrical  periodic  struc¬ 
tures  of  practical  dimensions  surrounding  the  area.  Po¬ 
tential  applications  include  planar  waveguide  devices  as 
well  as  cylindrically  symmetric  lasers  and  amplifiers  that 
operate  along  the  cylinder  axis. 

It  is  found  that  one  can  increase  the  region  of  source 
placement  for  which  significant  radiative  effects  occur  by 


pushing  the  Bragg  reflection  from  the  periodic  structure 
out  to  larger  radial  extents.  This  pushing  can  be  accom¬ 
plished  either  by  an  increase  in  the  inner  radius  of  the 
beginning  of  the  periodicity  or  by  use  of  a  long  (radially) 
grating  with  weak  modulation,  so  that  the  reflection  is 
spread  out  radially.  The  maximum  radial  extent  of  the 
Bragg  cavity  in  a  practical  structure  ought  to  be  limited 
only  by  propagation  losses. 

Finally,  an  analysis  of  the  harmonically  oscillating  cur¬ 
rent  line  source  frequency  shifts  that  are  caused  by  the 
presence  of  the  cylindrical  Bragg  cavity  finds  the  shifts  to 
be  small.  By  comparing  these  shifts  in  the  variation  of 
Bragg  reflectivity  of  the  cavity  with  frequency,  we  find 
that  the  source  remains  exactly  tuned  to  the  Bragg  fre¬ 
quency  to  a  good  approximation,  even  in  cases  of  strong 
enhancement  or  inhibition  of  radiation. 
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